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I. INTRODUCTION 



The detection by the COBE DMR instrument of fluctuations in the temperature 
distribution of the Cosmic Microwave Background Radiation (CMBR) on large an- 
gular scales is certainly one of the most significant cosmological results since the 
detection of the CMBR itself. These fluctuations provide valuable information about 
the nature of primordial perturbations believed responsible for the origin of structure 
in the Universe. The horizon radius at the epoch of last scattering of the CMBR 
corresponds to angular scales of about 2° on the sky, which implies that fluctuations 
on scales probed by COBE were not predominantly affected by causal processes or 
the nature of the matter constituents of the Universe at the time of last scattering 
of the CMBR. Indeed, the large-scale (greater than 2°) fluctuations arise from the 
Sachs- Wolfe effect when photons are either red or blue shifted as they climb out of, 
or fall into, gravitational potential wells 0. It is most likely that the fluctuations 
in the CMBR are the result of processes that occurred very early in the history of 
the Universe, so they yield vital information concerning the physics that led to the 
primordial perturbations. 

There are currently two very attractive scenarios for the origin of the primordial 
fluctuations: quantum effects during inflation, and gravitational effects of defects 
resulting from cosmological phase transitions. Both scenarios involve physics beyond 
the standard model of particle physics, involving energies in the range 10 10 GeV < 
E < 10 19 GeV, an energy scale we will refer to loosely as the Grand Unified Theory 
(GUT) scale. A major difference in the predictions of the two scenarios concern the 
Gaussian nature of the fluctuation pattern, and we should be able to use this to 
differentiate between the two possibilities in the near future. In this paper we will 
assume that the fluctuations are the result of inflation, and we discuss what might be 
learned about particle physics at very high energies from astronomical observations 
from which we can infer the primordial fluctuation spectrum. 

All models of inflation involve a period of rapid growth of the size of the Universe. 
This is most easily illustrated by considering a homogeneous, isotropic Universe with 
a flat Friedmann-Robert son- Walker (FRW) metric described in term of a scale factor 
a(t). Here, "rapid growth" means a positive value of a/a = — (4^(7^/3) (p + 3p) where 
p is the energy density and p the pressure. In all successful models of inflation, the 
Universe is dominated by some sort of scalar "potential" energy density V > that 
is positive, resulting in an effective equation of state p ~ —p ~ V, and hence a > 0. 
If one identifies the potential energy as arising from the potential of some scalar field 
4>, then (f) is known as the inflaton field. 

Even within this traditional view of inflation, there are two major ways to imple- 
ment the scenario. One way involves a first-order phase transition. In this method, 
either in the original proposal of Guth || or the latest version called extended in- 
flation ||, the inflaton is trapped in a meta-stable, or false- vacuum, state while the 
Universe inflates. Inflation is ended when the Universe undergoes a first-order phase 
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transition in which the inflaton field tunnels to its true-vacuum state. In the second 
method, inflation occurs because for some reason the inflaton field is displaced from 
its minimum and its potential energy density dominates the Universe; inflation occurs 
while the inflaton field is slowly evolving, or rolling, to its minimum || |j. It is this 
second class of "slow-roll" models we will consider in this paper.Q 

Although the early slow-roll models had potentials that were reasonably simple 
(Coleman- Weinberg, A0 4 , etc.), or at least polynomials in some scalar field, many 
attractive models have been developed where the scalar potential driving inflation 
is quite complicated. Perhaps the study of the density perturbations produced by 
inflation can shed some light on the nature of the potential. 

Broadly speaking, inflation predicts a very nearly Gaussian spectrum of density 
perturbations that is scale dependent, i.e., the amplitude of the perturbation depends 
upon the length scale. Such a dependence typically arises because the Hubble ex- 
pansion rate during the inflationary epoch in fact changes, albeit slowly, as the field 
driving the expansion rolls towards the minimum of the scalar potential. This implies 
that the amplitude of the fluctuations as they cross the Hubble radius will be weakly 
time-dependent. 

Within the context of slow-roll inflation, Hodges and Blumenthal have shown 
that any scale dependence for density perturbations is possible if one considers an 
arbitrary functional form for the inflaton potential, V(<j>). In this sense, inflation 
makes no unique prediction concerning the form of the spectrum and one is left with 
two options. Either one can aim to find a deeper physical principle that uniquely 
determines the potential, or observations that depend on V(<f)) can be employed to 
limit the number of possibilities. 

Improved observations of large-scale structure, of which COBE provides the most 
dramatic example at present, are important because they allow us, in principle, to 
determine the spectrum of primordial density perturbations. This may very well 
provide a direct experimental window on the physics of the Grand Unified era cor- 
responding to energy scales of the order 10 16 GeV. The purpose of the present work 
is to investigate to what extent information from the CMBR and large-scale galactic 
structure will allow us to reconstruct GUT physics. 

In the following section we will review the salient aspects of slow-roll inflation. 
In Section III we discuss the reconstruction of the inflaton potential from knowledge 
of scalar or tensor perturbations. Section IV illustrates the formalism by several 
examples in which the functional form of the potential is found from knowledge of 
the tensor and scalar perturbation spectra. Section V illustrates what can be learned 
about the potential from observations of the properties of the tensor and scalar spectra 
at a particular length scale. In Section VI the reader may find a discussion of how one 
determines the primordial density spectrum. Finally, Section VII offers an assessment 
of the prospectus for reconstruction of the inflaton potential. 

1 In reality, often the distinction between the two methods is not so clean, and it is possible to 
consider some types of first-order inflation models as variants of slow- roll models. See Ref. [[|. 
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II. REVIEW OF SLOW-ROLL INFLATION 



For the benefit of those not familiar with the generation of scalar and tensor 
perturbations in slow-roll inflation, we review the salient features in this section. 
Those comfortable with the basic results may wish to skip this section, and refer 
back to it as needed to understand notation and conventions. We set c = ft — 1, and 
define k 2 = 87iGn = 8n/m 2 Pl . 

Slow-roll inflation requires a scalar field to be displaced from the minimum of its 
potential at some time early in the evolution of the Universe. If during the evolution 
of the field to its minimum a region of the Universe is dominated by the potential 
energy of the field, then the volume of that region will undergo rapid expansion, 
inflate, and grow to encompass a volume large enough to contain all of the presently 
observed Universe. Eventually the potential energy ceases to dominate when the field 
evolves through a steep region of the potential and the field evolves so rapidly that 
the kinetic energy of the field comes to dominate. This is the end of inflation, and is 
followed by the scalar field oscillating about the minimum of its potential, with the 
inflaton field decaying and 're-heating' the Universe by conversion of vacuum energy 
to radiation. 

We are interested in the perturbations resulting from inflation. The "density" 
perturbations are usually described in term of fluctuations in the local value of the 
mass density. In a Universe with density field p(x) and mean mass density p , the 
density contrast is defined as 

*(x) = 6 -M = pW-P°. (2.1) 
Po Po 

It is convenient to express this contrast in terms of a Fourier expansion: 

<?(x) = A J^ exp(-ik • x)d 3 A;, (2.2) 

where A is simply some overall normalization constant, interesting only for those 
who enjoy keeping track of factors of 2tt. What is usually meant by the density 
perturbation on a scale A, (5p/p)\, is related to the square of the Fourier coefficients 
5 k : 



'8pY _ A ,k 3 \5 k \ 



p) x 



(2.3) 



where again we have included an overall normalization constant A'. The perturbations 
are normally taken to be (statistically) isotropic, in the sense that the expectation 
of 1 5k | 2 averaged over a large number of independent regions can depend only on 
k = |k|. The dependence of 5p/p as a function of A is the spectrum of the density 
perturbations. Of crucial importance is the relative size of scale A to the scale of the 
Hubble radius. The physical length between two points of coordinate separation d is 
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X(t) = a(t)d. A length scale comoving with the expansion will grow proportional to 
a(t). If a(t) < 0, as in the standard non-inflationary phase, then a(t) will grow slower 
than t. If a(t) > 0, as in the inflationary phase, then a(t) will grow faster than t. 

For a spatially flat isotropic Universe the Hubble expansion rate, H(t) = a/ a, is 
given by 3H 2 (t) = K 2 p(t). The inverse of the Hubble expansion rate, the Hubble 
radius A#(£) = H" 1 ^), is the scale beyond which causal processes no longer operate. 
In the non-inflationary phase A# increases linearly with time. Since in the non- 
inflationary phase Xh oc t, while X(t) increases more slowly than t, the Hubble radius 
increases faster than A(i), and a length scale A(i) will start larger than the Hubble 
radius (A > A#), cross the Hubble radius (A = A#), and then will remain inside the 
Hubble radius (A < A#). 

The story is different if we imagine that the Universe was once in an inflationary 
phase. In inflation H is roughly constant, so the Hubble length is roughly time inde- 
pendent. Thus, a given scale can start sub-Hubble radius, A < A#, then pass outside 
the Hubble radius during inflation, and then re-enter the Hubble radius after inflation. 
Thus, perturbations can be imparted on a given length scale in the inflationary era 
as that scale leaves the Hubble radius, and will be present as that scale re-enters the 
Hubble radius after inflation in the radiation-dominated or matter-dominated eras. 

Microphysics cannot affect the perturbation while it is outside the Hubble ra- 
dius, and the evolution of its amplitude is kinematical, unaffected by dissipation, the 
equation of state, instabilities, and the like. However, for super-Hubble radius sized 
perturbations one must take into account the freedom in the choice of the background 
reference space-time, i.e., the gauge ambiguities. As usual when confronted with such 
a problem, it is convenient to calculate a gauge-invariant quantity. For inflation it 
is convenient to study a quantity conventionally denoted ( ||. In the uniform Hub- 
ble constant gauge, at Hubble radius crossing ( is particularly simple, related to the 
background energy density and pressure p an d Po, and the perturbed energy density 
Pi- 

C = 5p/(p + Po ), (2.4) 

where 5p = p\ — po is the density perturbation. 

In the standard matter-dominated (MD) or radiation-dominated (RD) phase, ( 
at Hubble radius crossing (up to a factor of order unity) is equal to Sp/p. Thus, the 
amplitude of a density perturbation when it crosses back inside the Hubble radius 
after inflation, (5p/p)noK^ is given by ( at the time the fluctuation crossed outside 
the Hubble radius during inflation. 

As inferred from the adoption of (, the convenient specification of the amplitude of 
density perturbations on a particular scale is when that particular scale just enters the 
Hubble radius, denoted as (5p/p)hor- Specifying the amplitude of the perturbation 
at Hubble radius crossing evades the subtleties associated with the gauge freedom, 

2 The notation "HOR" follows because often in the literature the Hubble radius is referred to 
(incorrectly) as the horizon. 



4 



and has the simple Newtonian interpretation as the amplitude of the perturbation in 
the gravitational potential. Of course, when one specifies the fluctuation spectrum at 
Hubble radius crossing, the amplitudes for different lengths are specified at different 
times. 

Now let us turn to the scalar field dynamics during inflation. Consider a minimally 
coupled, spatially homogeneous scalar field 0, with Lagrangian density 

c = - v{4>) = <j> 2 /2 - v ((/>). (2.5) 

With the assumption that is spatially homogeneous, the stress-energy tensor takes 
the form of a perfect fluid, with energy density and pressure given by = <fi 2 /2+V((f)), 
and = 4> 2 /2 — V(<f>). The classical equation of motion for is 

+ 3#0 + K'(0) = 0, (2.6) 

with the expansion rate in a flat FRW spacetime given by 

H 2 = y Q0 2 + V{4>)) . (2.7) 

Here dot and prime denote differentiation with respect to cosmic time and respec- 
tively. We assume that inflation has already provided us with a flat universe by the 
time the largest observable scales cross the Hubble radius. 

By differentiating Eq. fl2.7| ) with respect to t and substituting in Eq. (|2.6|), we 
arrive at the "momentum" equation 

2H = -k 2 2 . (2.8) 

All minimal slow-roll models are examples of sub-inflationary behavior, which is de- 
fined by the condition H < 0. Super-inflation, where H > 0, cannot occur here, 



though it is possible in more complex scenarios || [10| . We may divide both sides 
of this equation by if this quantity does not pass through zero. This allows us to 
eliminate the time-dependence in the Friedmann equation [Eq. (|2.7f) l and derive the 
first-order, non-linear differential equations 

(H'f-l^H 2 = -^V(0) (2.9) 
k 2 = -2H'. (2.10) 

A common framework for discussion of inflation is the slow-roll approximation, 
though let us emphasize here that in much of our treatment of inflaton dynamics we 
shall not need to resort to it. We can define two parameters, which we will denote as 
slow-roll parameters, by 3 




3 These definitions differ slightly from, and indeed improve upon, those of Refs. ([il|), ( ^2[ ) which 
were made using the potential rather than the Hubble parameter. As defined here they possess 
rather more elegant properties. 
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Slow-roll corresponds to {e, \r}\} <C 1. These conditions correspond respectively to 
the cases when the first term in Eq. (|2.9|) and the first term in its 0-derivative can be 
neglected. 

With these definitions, the end of inflation is given exactly by e = 1. A small 
value of r] guarantees 

3H(t)~ -V ((/)), (2.12) 

which is often called the slow-roll equation.^ Although the terminology "slow-roll 
approximation" is normally used rather loosely, one could imagine carrying out a 
formalized perturbation expansion in the slow-roll parameters, and we shall refer to 
such results later. 

Density perturbations arise as the result of quantum-mechanical fluctuations of 
fields in de Sitter space. First, let's consider scalar density fluctuations. To a good 
approximation we may treat the inflaton field as a massless, minimally coupled 
field. (Of course the inflaton does have a mass, but inflation operates when the 
field is evolving through a "flat" region of the potential.) Just as fluctuations in the 
density field may be expanded in a Fourier series as in Eq. fl2TT|), the fluctuations in 
the inflaton field may be expanded in terms of its Fourier coefficients £0k : <50(x) oc 
/ 50k exp(— ik ■ x)d 3 k. During inflation there is an event horizon as in de Sitter space, 
and quantum-mechanical fluctuations in the Fourier components of the inflaton field 
are given by [[HJ 

fc 3 |50 k | 2 /27r 2 = (#/27r) 2 , (2.13) 

where H/2tt plays a role similar to the Hawking temperature of black holes. Thus, 
when a given mode of the inflaton field leaves the Hubble radius during inflation, it 
has impressed upon it quantum mechanical fluctuations. In analogy to Eq. (|2.3|) , what 
is called the fluctuations in the inflaton field on scale k is proportional to fc 3//2 |<50k|, 
which by Eq. ( |2.13|) is proportional to H/2n. Fluctuations in lead to perturbations 
in the energy density: 

5 P4 , = 5<P{dV/d<P). (2.14) 

Now considering the fluctuations as a particular mode leaves the Hubble radius during 
inflation, we may construct the gauge invariant quantity ( from Eq. ( |2.4[ ) using the 
fact that during inflation p + p = 2 : 



4 With the definition of e in Refs. (p"i"|), (|l2|), this result is true only in the slow-roll approximation. 

5 Note the difference between slow-roll inflation and the slow-roll equation. Slow-roll inflation is 
a model where inflation occurs where the scalar field is slowly evolving to its minimum, while the 
slow-roll equation implies that <p can be neglected. 
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Now using Eq. ( |2.9|) and Eq. ( |2.10|) , the amplitude of the density perturbation 
when it crosses the Hubble radius after inflation is 

if), ^T2 AM ^ww\ x ^ywr (2 - 16) 

where H(4>) and H'{4>) are to be evaluated when the scale A crossed the Hubble 
radius during inflation. The constant m equals 2/5 or 4 if the perturbation re-enters 
during the matter or radiation dominated eras respectively.^ Now we wish to know 
the A-dependence of (5p/p)\, while the right-hand side of the equation is a function 
of cf) when A crossed the Hubble radius during inflation. We may find the value of the 
scalar field when the scale A goes outside the Hubble radius in terms of the number 
of e-foldings of growth in the scale factor between Hubble radius crossing and the end 
of inflation. 

It is quite a simple matter to calculate the number of e-foldings of growth in the 
scale factor that occur as the scalar field rolls from a particular value <fi to the end of 
inflation e : 

w-tZ (2 - l7) 

The slow-roll conditions guarantee a large number of e-foldings. The total amount of 
inflation is given by iVroT = N((j)i), where fa is the initial value of <fi at the start of 
inflation (when a first becomes positive). In general, the number of e- folds between 
when a length scale A crossed the Hubble radius during inflation and the end of 
inflation is given by |jTT] 



N(X) = 45 + ln(A/Mpc) + - ln(M/10 14 GeV) + - ln(T RH /10 10 GeV), (2.18) 

3 3 

where M is the mass scale associated with the potential and Trh is the "re-heat" tem- 
perature. Relating N(\) and N(4>) from Eq. Q2.17Q results in an expression between 
(ft and A. Hopefully this dry formalism will become clear in the example discussed 
below. 

In addition to the scalar density perturbations caused by de Sitter fluctuations in 
the inflaton field, there are gravitational mode perturbations, g^ u — > <?™ w + h^, 
caused by de Sitter fluctuations in the metric tensor JTJjrrSJ. Here, ff^ w is the 
Friedmann-Robertson-Walker metric and are the metric perturbations. That 
de Sitter space fluctuations should lead to fluctuations in the metric tensor is not 



6 The 4 for radiation is appropriate to the uniform Hubble constant gauge. One occasionally 
sees a value 4/9 instead which is appropriate to the synchronous gauge. The matter domination 
factor is the same in either case. Note also that it is exact for matter domination, but for radiation 
domination it is only strictly true for modes much larger than the Hubble radius, and there will be 
corrections in the extrapolation down to the size of the Hubble radius. 
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surprising, since after all, gravitons are the propagating modes associated with trans- 
verse, traceless metric perturbations, and they too behave as minimally coupled scalar 
fields. The dimensionless tensor metric perturbations can be expressed in terms of 
two graviton modes we will denote as h. Performing a Fourier decomposition of h, 
/i(x) oc / <5/ifcexp(— ik • 5i)d 3 k, we can use the formalism for scalar field perturbations 
simply by the identification <50 k — * h^/ny/2, with resulting quantum fluctuations [cf. 
Eq. 

k 3 \h k \ 2 /2ir 2 = 2k 2 (H/2ti) 2 . (2.19) 

While outside the Hubble radius, the amplitude of a given mode remains constant, 
so the amplitude of the dimensionless strain on scale A when it crosses the Hubble 
radius after inflation is given by 



H ° R ^ O (0) = ^^)~^, (2.20) 



m 



2 

PI 



where once again H(4>) is to be evaluated when the scale A crossed the Hubble radius 
during inflation. 

As usual, it is convenient to illustrate the general features of inflation in the 
context of the simplest model, chaotic inflation ||, which is to inflationary cosmology 
what drosophila is to genetics. In chaotic inflation the inflaton potential is usually 
taken to have a simple polynomial form such as V(4>) = A0 4 , or V(<p) = fJ?4> 2 - 
For a concrete example, let us consider the simplest chaotic inflation model, with 
potential V{4>) = /J, 2 (f) 2 pf. This model can be adequately solved in the slow-roll 
approximation, yielding 

<f>{t) = 




a(t) = di exp — p I (pit ^— t 2 



y/3. 



K 



with inflation ending at K(p e = \[2 as determined by e = 1, where e was defined in 
Eq. ( |2.11| ). The number of e-foldings between a scalar field value 0, and the end of 



inflation is just 

c 2 f V2/K K 2^2 l 



w)=--r &^=— - - ( 2 - 22 ) 

yv! 2 J4> H'ty) ^42 v 1 

Equating Eq. ( |2.22| ) and Eq. ( |2.18| ) relates and A in this model for inflation: 

K 2 (j) 2 /4 = [45.5 + ln(A/Mpc)] . (2.23) 
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Using Eq. ( |2.16|) and Eq. ( p.20|) , As and A G are found to be 



A S (X) = (V2k^/VV2^) [45.5 + ln(A/Mpc)] 

A G {\) = (K/i/v / 127 r ) [45.5 + ln(A/Mpc)] 1/2 . (2.24) 

We can note three features that are common to a large number of (but not all) 
inflationary models. First, As and A G have different functional dependences upon A. 
Second, A G and A s increase with A. Finally, A s > A G , for scales of interest, although 
not by an enormous factor. 

To conclude this exercise, it is worth reminding the reader how little of the inflaton 
potential is available for reconstruction. The scales of cosmological interest at the 
present epoch lie in the range l/i -1 Mpc for galaxies up to the current horizon size 
of 6000/i _1 Mpc, where as usual h is Hubble's constant in units of 100 km s _1 Mpc _1 . 
Taking the present horizon distance to have crossed the Hubble radius 60 e-foldings 
from the end of inflation, we see that we only sample the small region of the potential 
V(4>) for G [2. 7m pi, 3. 0m pi]. By any standards, the density perturbations from 
inflation we can actually sample represent an extremely small region of the potential. 
However it should be realized that although we have potential information about a 
small region of the potential, any information about the GUT potential, no matter 
how meager, is precious! Indeed, in the exploration of GUTs, cosmology may reveal 
the first "piece of the action." 

We have one further piece of information, which is that we know that the remain- 
der of inflation must occur in the remaining section of the potential, with the scalar 
field coming to rest with V(4>) = 0. Although this represents a significant constraint 
on the potential on scales below those that large-scale structure observations can sam- 
ple, it does still leave an uncountable infinity of possible forms in this region. (One 
other constraint in this region comes from primordial black holes, whose abundance 
can in principle be calculated from the spectrum. Should black hole formation be 
copious, this constrains the spectrum at the mass scales corresponding to the size of 
detectable black holes, which are most prominent at around 10 15 g |[L7|| .) 

To conclude this section we call the reader's attention to Figure 1, which illustrates 
the procedure. The figure illustrates a scalar field <fi rolling down a potential V((f>). 
At some point in the evolution the slow-roll conditions break down and inflation 
ends. We can count back from this point the number of e-foldings from the end of 
inflation, and use this information to find a relationship between the value of <fi in 
the evolution, and the length scale A leaving the Hubble radius at that point. While 
evolves, quantum fluctuations imprint scalar and gravity-wave perturbations upon 
each scale as it leaves the Hubble radius. The scalar perturbations depend upon 
the potential and its derivative, while the gravitational modes depend only upon the 
potential. In principle, A G and As are probed by observations of large-scale structure 
and by measurements of CMBR fluctuations. The length scale, and corresponding 
angular scales, of several important observations are indicated. 
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In the next section we will discuss the procedure for reversing the process discussed 
in this section; knowing A G (X) and A S (X), how does one determine V(0)? 



III. RECONSTRUCTION OF THE POTENTIAL 



A number of authors fll2", n^l9| , p0| | have recently emphasized the possibility that 



tensor modes excited during inflation, corresponding to gravitational waves, may play 
an important role in generating microwave background anisotropies. We thus develop 
an extension of the potential reconstruction methods of Hodges and Blumenthal |Q 
to include tensor as well as scalar modes. As discussed in the previous section, the 
expressions for the amplitudes of the scalar and tensor modes may be written as 

M4>) - ^ W 



8vr 3 / 2 |if'(0)| 



m*) = (3.1; 



respectively. Note that the definition of As ((f)) in Eq. (|3.1| ) is related to the power 
spectrum P 1//2 (/c) defined in Hodges and Blumenthal |7j by 



P x ' 2 (k) = 3V2^A s ((f)). (3.2) 

Utilizing the slow-roll approximation, there are useful expressions for the scale- 
dependence of the spectra, the spectral indices, to first-order in departure from slow- 
roll. These are 

l-n= d\n[A s (X)]/d\nX = 4e* - 2^ 

n G = d]n[A 2 G (X))/dlnX = 2e* (3.3) 

where "*" indicates evaluation at the time when the scale A passes outside the Hubble 
radius during inflation. In keeping with convention we drop the subscript S on the 
scalar mode index. Whenever the slow-roll conditions are closely obeyed, the spec- 
trum is close to scale invariant. When this is not true, there are in general corrections 
to the expressions for the fluctuations at the next order in an expansion in slow-roll 
parameters. 

The reader may have noticed that although we are keeping the equations of motion 
general (i.e., not subject to a slow-roll approximation), our expression for the scalar 
modes in Eq. (|3.1|) is an expression based on the slow-roll approximation, {e, \r]\} <C 1. 
Ideally, one would like to completely abandon the slow-roll regime, because within it, 
the scalar spectrum is always close to the scale-invariant case and the gravitational 
wave amplitude is always small, as we have seen. In practice, it seems very possible 
that should inflation have occurred, it may well have been pushing the outside of the 
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slow-roll approximation envelope, and indeed much of the recent interest has been 
in the possibilities of both tilt and gravitational waves. True reconstruction assumes 
nothing about V(4>) (flatness, etc.) except that it inflates. Unfortunately, although 
we are able to keep the dynamics completely general, general expressions are not 
available for the perturbation spectra. 

Recently, an improvement has become available in the form of general expressions 
for the spectra to first-order in departure from slow- roll [2~I|. These give rise to 
"first-order corrected" spectra, which can be written 



I uncorr 



AT = [l-e + (2-ln2- 7 )(2e-77)]^ 
A C ° TT = [1 + (1 -In 2 -7)e]A™ r , (3.4) 

where 7 ~ 0.577 is Euler's constant. If slow-roll is breaking then these can rep- 
resent a significant improvement on the uncorrected results, but unfortunately the 
reconstruction loses its analytic tractability. The one exception to this is the case of 
power-law inflation — in that case the effects of the corrections cancel exactly |22|] in 
the reconstruction equation Eq. ( |3.10 ) we derive below. 



(3.5) 



Rather than resort immediately to numerical construction, we elect instead to 
make the operational choice that we shall adopt the slow-roll expressions for the 
spectra. A reconstruction can then be made subject to a consistency check that the 
slow-roll conditions are indeed satisfied; if not, then our formalism will have to be 
enhanced to incorporate these improvements. 

It is clear that the ratio of amplitudes of the scalar and tensor modes is given by 

A G s/2 \H'\ r 

and if e <C 1, then Aq/As < 1. It is possible that the COBE satellite is in fact 
observing a sum of contributions from the tensor and scalar fluctuations, as opposed 
to the pure scalar modes as originally thought. If these are uncorrelated and obey 
Gaussian statistics, the quantity of observational interest on large angular scales is 
the sum of the squares[] 
2 

Tfl 

S\<P) = —AU<P) + AU<P). (3.6) 

Using m = 2/5 and recalling that e must be less than unity, we see immediately 
that the tensor modes dominate S 2 (<f)) if 2/25 < e < 1, or equivalently if k/5 < 
\H'\/H < k/s/2. The largest relative tensor contribution to S 2 ((f)) obtains for e = 1: 
2A 2 j/m 2 A 2 s = 25/2 for m = 2/5.0 Although it is not mandatory (one can break 

7 The relative weighting of A s and A 2 G in this equation is that appropriate to large angle 
anisotropies (greater than 2°) in the slow-roll approximation. This is discussed in depth in Sec- 
tion VI, and exact weighting formulae provided there. 

8 If one is performing a theoretical reconstruction of the potential by specifying either Aq or 
Ag, it is essential to ensure this condition is always satisfied for consistency. Indeed, observations 
violating Aq/As < 1 would immediately rule out the models we are considering. 
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slow-roll only in the 77 parameter as in natural inflation), the gravitational wave 
contribution is typically significant whenever there is a deviation from the slow-roll 
regime. 

On the other hand, the gravitational waves behave as relativistic matter when they 
re-enter the Hubble radius and do not interact with the other matter components. 
Consequently their energy density redshifts as a -4 , which implies that only scalar 
modes affect the CMBR anisotropy on angular scales 6^2°. However the anisotropy 
on these scales is also affected by the form of dark matter present. [For a recent 
discussion of some of these issues, see Ref. (J2"5|).] 

To proceed, we shall assume that the functional forms of As(X) and Aq(X) are 
known explicitly and defer until Section VI a discussion on the many difficulties 
associated with determining these quantities from observation. Our initial aim is 
to develop a framework which allows the inflaton potential to be determined. We 
consider general inflationary behavior for the field equations ( |2.9| ) and ( |2.10| ) and it 
proves convenient to parameterize the full set of solutions in terms of the function 
if (A(0)), where A is the scale. Eqs. ( |3.1| ) now become 

V2k 2 



A S (X) 
A G (X) 



^3/2 

K 



H\X) 
H(X). 



dX d(j) 
dHdX 



(3.7) 



47r 3/2 

Each length scale A is associated with a unique value of <p when that scale crossed 
the Hubble radius during inflation. We will indicate that relationship by writing A(0). 
Now when a present length scale A crossed the Horizon radius during inflation with 
scalar field value 0, its physical size was if _1 (0). The physical size grew between 
horizon crossing and today, and is now simply A(0) = i/ _1 (0)ao/a(0), where ao 
is the present value of the scale factor and a(<p) was the value of the scale factor 
when the scale crossed the Hubble radius during inflation. Now we can make use of 
Eq. ( |2.17| ) to relate a(<f)) to the value of the scale factor at the end of inflation, a e : 
a(<p) = a e exp[— N(<f>)]. This allows us to express X(<p) as 

exp[N(4>)] a 



A(0) 



H{4>) 



where N(<f>) is given by Eq. ( |2.17 ). Differentiating Eq. (|3.8| ) with respect to 



dx(<f>) _ ± 

d<p 


k_ (As_ _ 


A G \ 




As) 


and taking the ratio of Eqs. 




K_ Ac _ 


dX dlxiAc 




V2 A s 


dip dX 





A. 



(3.8) 
yields 
(3.9) 

(3.10) 



9 If one were to use the 'first-order corrected' expressions for the spectra discussed earlier, the 
right hand side of Eq. (|3.10|) would be multiplied by (1 - 1.27e + 1.27t?). 
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Note that expression (2.10) in Hodges and Blumenthal JF| consists of only our first 
term in Eq. (|3.9| ), indicating their assumption of slow- roll behavior. Substituting Eq. 
( PI) into Eq. ( CT ) gives 

A dA G (A) A 2 G (A) 



A G (A) rfA ,4f(A)-A G (A) 



(3.H) 



Note that the left hand side is just equal to n G /2. This equation is similar to Eq. (9) 
in Davis et al. [[HJ , provided one interprets their n as being the tensor index and not 
the scalar one. It clearly shows that there exists a correspondence between the scalar 
and tensor modes and is valid for an arbitrary interaction potential. In principle, if 
the scale dependence of either the scalar or tensor modes is known, the other can be 
determined from Eq. ( |3.11| ). If only A G (X) is known, then As(X) follows immediately 
by differentiation. However, if only As(X) is known, a first-order differential equation 
must be solved to find the form of Aq(X). Thus, knowledge of only the scalar spectrum 
leaves an undetermined constant in the tensor spectrum. 

Once the form of the tensor spectrum is known, the potential, as parametrized by 
A, may be derived by substituting Eqs. ([3.7|) into Eq. ( |3.10| ). We find 



v[4>(\)) 



16tt 3 A g (A) 



« 4 



(3.12) 



Am. 

Finally, integration of Eq. ( p.9| ) yields the function = 0(A) given by 
j.f \ \ | ^ f x dX' A S (X')A G (X') 

m = ± ~J -y au\') - Auyy (3 - 13) 

We have absorbed the integration constant by taking advantage of the freedom to 
shift by a constant. The functional form of V(0) follows by inverting Eq. ( [3.13 ) 



and substituting the result into Eq. ( |3.12[ ). It will also prove convenient at times to 



express in terms of A G . If the functional form of As as a function of A G is known, 
A S [A G ], then using Eq. ( gTj ) in Eq. (|3TT3| ) gives 

i^r^WU. ,3,4) 



K 



The reconstruction equations are Eqs. ( p.ll|) , (|3.12|) , and (|3.13| ). It is worth 



emphasizing again that for any choice of A G (X), there is a unique associated As(X) 
and V{4>) (at least in the slow-roll approximation), but that the converse is not true. 
As shown by Hodges and Blumenthal |7J, the scalar spectrum leaves an undetermined 
constant in the tensor spectrum, and as the equation relating V and Aq is non-linear, 
different choices of this constant might lead to functionally different forms of the 
potential [pl |. In order to reconstruct the potential from scalar modes, one needs 



an additional piece of information. Technically what is needed is knowledge of the 
functional dependence of As upon A G , As[Aq]. This can be fixed either by knowledge 
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of the amplitude of the tensor spectrum at a single scale, which would fix A G uniquely, 
or knowledge that A s is independent of A G . As A G cannot be independent of A, the 
latter possibility arises only if As(A) is constant. 

It is also worth emphasising consistency, which can provide an important check. If 
our inflationary assumptions are correct, then the two spectra are intimately related 
as illustrated above. However, observations are typically subject to both systematic 
and statistical errors, and within these one might find that measured spectra are not 
exactly consistent. Were one to be confronted with such data, one would like some 
prescription by which to decide how to best reconcile the data, to generate some kind 
of "maximum likelihood" reconstruction. Such a procedure would presumably also 
allow one to demonstrate that the measured spectra were not compatible with each 
other within the inflationary paradigm, if indeed inflation were not the correct source 
of the fluctuations. In practice, the situation is skewed by the scalar fluctuations being 
considerably easier to observe than their tensor counterparts, and it seems prudent 
to await the arrival of considerably better data before properly contemplating how 
one would deal with the possibility of marginally incompatible observations. 

The reconstruction procedure simplifies if Aq(\) <C A s (X) (i.e., e <C 1): 

A dA G (X) A%{X) 



A G (X) dX A|(A) 

V2 f x d\' A G (X') 



0(A) = ±- 



k J X' A S (X') 



V[<K\)\ = ^?^g(A). (3.15) 



K 4 



We conclude this section by summarizing the conditions necessary for the pertur- 
bation amplitudes to increase or decrease with increasing wavelength. Such informa- 
tion alone can place strong limits on the functional form of the potential. The scales 
that first cross the Hubble radius are the last to re-enter during the radiation or mat- 
ter dominated eras (see Fig. 1). Consequently, the amplitudes of the modes increase 
(decrease) with wavelength if they decrease (increase) with time during inflation. 
Immediately we conclude that 

dA G 

it > < 316 > 

for all sub-inflationary (H < 0) models. One requires an era of super-inflation (H > 0) 
if this inequality is to be reversed. Super-inflation is only possible with a minimally 
coupled self-interacting scalar field if the spatial hypersurfaces of the manifold have 
positive-definite curvature 0. An observation indicating dA G /dX < would there- 
fore require some of the main assumptions made in the inflationary analysis to be 
significantly altered. Within the context of the FRW Universe, for example, one 
would need to extend the gravitational sector of the theory beyond general relativ- 
ity, or assume that the value of the density parameter was significantly larger than 
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unity at first Hubble radius crossing. Indeed, Eq. (|3.16|) implies that any effects of 
the gravitational waves on the CMBR anisotropy will always be enhanced on larger 
angular scales in the models considered here. 

On the other hand, it is possible for the scalar spectrum to decrease with wave- 
length. By writing dAs/dX = {dAs / d(j)){d(j) / dX) and employing Eqs. ( |3.7p and ( |3.9| ), 
one finds that a necessary and sufficient condition for scalar modes to be decreasing 
in amplitude with increasing wavelength is 

(3.17) 

In terms of the slow-roll parameters, this can be written as 2e < rj. As e is positive 
by definition, this condition is not easy to satisfy, particularly in the late stages 
of inflation where e must increase towards unity. A necessary, but not sufficient, 
condition for Eq. (|3.17| ) to hold is that the potential be convex, V" > 0. Therefore, if 




the field is located near a local maximum of the potential, as in natural inflation [25 
for example, the amplitude will always increase with A. 

In conclusion, it is clear that any scale dependence for the spectrum of gravita- 
tional waves is possible in principle, subject to condition ( |3.16| ). Secondly, the most 



useful parameter mathematically in the reconstruction process is A G (X), because once 
this is known the potential can be derived in a rather straightforward manner. 



IV. RECONSTRUCTING THE FULL POTENTIAL 

Before proceeding to analyze the possibilities for obtaining the spectra observa- 
tionally, we shall first illustrate some examples in reconstruction in order to demon- 
strate the techniques. We shall examine four cases of increasing complexity. These 
four cases will reconstruct to familiar potentials. 

A. Polynomial potentials 

Let us first reconstruct the n 2 <ft 2 chaotic potential model worked out in Section II. 
We will then generalize the result for construction of polynomial potentials. 

Recall that using the slow-roll approximation for the potential V(4>) = fi 2 <p 2 we 
found perturbation spectra A G (X) = ct[(3 + ln(A/Ao)] 1//2 and As{\) = \/2A G (\)/a, 
with a 2 = K 2 /x 2 /127r 3 , f3 = 45.5, and Ao = 1 Mpc. We must keep in mind that 
these solutions were obtained in the slow-roll approximation. Since the slow-roll 
approximation implies that A G <C As, we must reconstruct using Eqs. (|3.15|) . 

First, let us reconstruct assuming that observations provide two pieces of infor- 
mation: A G (\) is of the form A G (X) = a[/3 + m(A/A )] 1/2 , and A G (X) < As (A). Then 
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the differential equation for oL4g(A) / 'd\ in Eq. ( |3.15|) can be used to yield a unique 
scalar spectrum, As(A) = y/2A G (X) / a, as anticipated from the calculation in Section 
II. (Of course As(A) could be found without the assumption that A G (\) <C As (A), 
but it would be different.) 

Now that we know both A?(A) and As(A), we can find 0(A) from the second 
equation in Eq. (|3.15|) : 

(3 + ln(A/A ) = k 2 <P 2 /A. (4.1) 

Finally, we can use the last equation in the slow-roll reconstruction procedure to 
give 

V(<t>) = ^A%(X) = ^a 2 [P + ln(A/A )] = ^a 2 2 . (4.2) 

Exactly as expected, the potential is of the form V(<p) = fi 2 (p 2 , with /i 2 = YItv'q 2 / k 2 . 
Thus, we have successfully reconstructed the potential. 

We began with the assumption that A G (X) is known. If we had started with the 
assumption that the scalar spectrum is known and of the form As (A) = y/2a[/3 + 
ln(A/Ao)], the differential equation for A?(A) would give 

A- G 2 (\) = a~ 2 [P + IntA/Ao)]- 1 + C, (4.3) 

where C is arbitrary. Fixing As*(Ao) = a/3 1 ^ 2 fixes C — 0, and reconstruction would 
proceed exactly as before. Other choices of C would lead to different potentials, with 
different predictions for Aq. 

Now let's consider a slightly more general tensor mode spectrum: Aq = a[(3 + 
ln(A/A )] 7 with 7 = constant, again with Aq{\) ^ As(A). The differential equation 
for dAc(\) / d\ gives 

A?(A) = [a/V7] IP + ln(A/A )] (27+1)/2 . (4.4) 



The solution for 0(A) is the same as Eq. (|4.1|) with k — > n/y/QFj. Using this in the 
reconstruction of the potential gives 

\r(j.\ 487:3 a2 { ^ 48tt 3 a 2 ^ 

= ^T-^c(A) = jj^4> 7 - (4-5) 

An oft studied case is 7 = 1, which reconstructs to V(</>) = A0 4 with scalar and tensor 
perturbations 

A G (X) = a[/? + ln(A/A )] 

A S (X) = «[/? + ln(A/A )] 3/2 . (4.6) 



B. Harrison— Zel'dovich potentials 
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Let us now look at potentials which give rise to the Harrison-Zel'dovich spectrum, 
A S (X) = as = constant. Such spectra are actually rather unlikely; most inflationary 
models exhibit a decrease in amplitude with decreasing scale which is significant now 
given the accuracy of observations. 

We start reconstruction by considering the differential equation relating A G and 
A s [Eq. ( full) : 



A dA G (\) A 2 G {\) 



A G {\) dX a 2 s -A 2 G {\) 
which has solution 



(4.7) 



4 



ln(A/A ) = -f (^^y - -grj - \n(A G (X)/A G0 ). (4.8) 

In general there is no closed- form expression for A G (X). 

We can reconstruct the potential in two steps. Since As is a constant, we can find 
A G in terms of <fi by Eq. ( |3.14| ): 

A 2 G {<f>) = 2a|/«V 2 - (4-9) 
Now we can substitute this into the equation for V in Eq. (|3.12 ) to give 



v{4>) 



,3 



167T 



2 



1 1 



(0/0) 2 u/w 



(4.10) 



where 4> — V%/ k>. 

It should be emphasized that this is the only inflaton potential which leads to an 
exactly scale-invariant spectrum of scalar density fluctuations. It arises as a special 
case of "intermediate" inflation where the scale factor expands as a oc exp(t^) 
with < / < 1; the above potential corresponds to choosing / = 2/3. In contrast, 
the spectrum of gravitational waves is not scale invariant. It is generally true that 
inflation cannot lead to scalar and tensor perturbation spectra that are both constant 
in A. It is interesting to note that potentials of this form arise when supersymmetry 



is spontaneously broken [27| 



We can reinterpret these results in terms of the slow-roll parameters. It is clear 
that to obtain a flat spectrum we require 2e = r), but e and 77 are not determined 
separately. There are some interesting limiting cases. If we allow as to tend to 
infinity, this corresponds to e tending to zero. In this limit the potential becomes 
flat, with its constant value being that which gives the desired gravitational wave 
amplitude. As as is reduced from infinity, e increases away from zero preserving 2e = 
77. Once e becomes big enough, there will be slow-roll corrections which destroy the 
flatness of the spectrum. It is interesting to note that although slow-roll automatically 
guarantees a spectrum which is close to flat, it is perfectly possible for a spectrum 
close to flatness to arise when the slow-roll conditions are not well obeyed. 
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Scalar 


Small gravitational 


Large gravitational 


Spectrum 


wave contribution 


wave contribution 


Nearly flat 


Polynomial 


Harrison- Zel 'dovich 


spectrum 


potentials 


potential 


Tilted 


Hyperbolic 


Exponential 


spectrum 


potentials 


potentials 



Table 1: Possible behaviors for spectra in several inflationary models. 



These potentials, which exhibit little tilt but which can have substantial gravita- 
tional waves, are also of interest in that they complete a square of possible behaviors 
in different inflationary models, as shown in Table 1. Indeed, such a model performs 
well on most large-scale structure data with the exception of intermediate-scale galaxy 
clustering data. 



C. Exponential potentials 



Generalizing away from the flat scalar spectrum, the simplest (and possibly most 
likely) case is where the amplitudes have a simple power-law dependence, 

A S (X) = a s (\/\ o y, u^O, (4.11) 

where as is a constant. The recent measurements from COBE [0 alone provide the 
constraint —0.3 < v < 0.2 at the 1-sigma level. Incorporating specific choices of dark 
matter and including clustering data allows one to do better; for instance in a cold 
dark matter model (CDM) it has been shown ]2S| that v < 0.15 at 95% confidence in 
models with no gravitational waves, and v < 0.08 (again 95% confidence) in power-law 
inflation which does have significant gravitational wave production. 
Aq{X) satisfies the differential equation Eq. Q3.11|) 



A dA G (X) _ Al{\) 



A G {\) dX a 2 s (X/X y--A 2 G (\y 

Obtaining the general form for Aq (A) is difficult. However there are some specific 
solutions which are of interest in that they relate to known examples of inflationary 
potentials. One obvious solution to Eq. ( f4.12p is A G (X) = a 9 (A/A ) iy , with 
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Note that in this simple case, Ac /As = a g /as, a constant independent of scale, but 
that as v — > the magnitude of the tensor contribution reduces significantly. We can 



simply integrate Eq. ( |3.13| ) to obtain 

(4.14) 




Substituting this expression into Eq. ( gjjp gives the final result 



7(0) = Voexp(±0/$, (4.15) 

with 



16n 3 a 2 g 2z/ + 3 - x I 2u 
v o = — —Ti <t> = K \—T- 416 

K V + 1 V V + 1 

Thus we see that a power-law behavior for the amplitude of the scalar and tensor 
modes is obtained from an exponential potential, and is therefore consistent with 
power-law models of inflation [K| . 



It is interesting to note that this result for cf> coincides with the exact result for 
power-law inflation, whereas if slow- roll were strictly applied one would get = 
Ky/2is, being the above to lowest order in v. Thus our hybrid of general equations of 
motion but slow-roll spectrum definitions certainly offers improved results over the 
usual slow-roll method in this case. 

Note that as v — > +oo the relative slope of the potential, as determined by a, 
becomes independent of v. The limit a{y = oo) = \/2k corresponds to the Milne 
Universe a(t) oc t and represents the limiting solution for inflation to occur. As v is 
increased in Eq. (|4.15 ) the only real effect is to increase the height of the potential 



through the Vq term. 

Rather than the equal power behavior, consider the more general example 

A S (X) = a s (\/\ o y ; A G (X) = a G (A/A o r, (4.17) 

where as and ac are constants. It is trivial to show that these spectra are solutions to 
Eq. (|3 . 1 1|) or the differential equation in Eq. ( p. 151 ) only if a = v. Thus, observation 



of spectra that are exact power-laws with different powers would rule out the class of 
inflationary models we consider as the source of the perturbations. 

D. Hyperbolic potentials 

Let us return to the differential equation for Aq{\) in Eq. ( |4.12j ), but in the 
AgW *C A s (X) limit. The equation becomes 



dA G = A% 1 ^ M ^_,„ 
dX X a 2 



(A/A )-^. (4.18) 

s 
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This equation has general solution in terms of an undetermined constant (3: 

(A/Ao) 2 ' 



/3(A/A 



\2u ' 



(4.19) 



We will see that different functional forms for the potential reconstruct depending 
upon the sign of (3. Of course f3 can be determined by measurement of Aq on any 
one scale. 

As (3 —>■ we recover the power-law spectra for Aq(X) and A$(X) with equal power- 
law slopes. This case was just considered above. For small scales, \(3\ (X/X ) 2u <C 1, 
we also recover the above case of equal power-law slopes for either choice of the sign 
of (3. For (3 > we can take the limit of large scales, j3{X/Xo) 2v ^> 1, in which case 
Ag(X) asymptotically approaches a constant. 

Recall that in the A G <C A s reconstruction procedure, Eq. ( |3.15| ) gives 



V[i>(X)) 



48tt 



(4.20) 



Now we must find 0(A) and invert. 

The integral expression for 0(A) from Eq. ( |3.15 ) is 



2;/ 



a dX! 



X' [1±|/3|(A'/A ) 2 1 1/2 ' 



(4.21) 



with "+" for positive (3 and "— " for negative j3. The constant is the same as that of 



the previous subsection (in the slow- roll approximation), 
to the integral are 



-2Arccsch 



-2Arcsech 



V^(A/A r P>0 



(X/XoT (3<0. 



K\l2v. The solutions 



(4.22) 



These expressions are easily inverted to give A(0), and the potential reconstructs to 
cosh" 2 (0/20) (3 > 



^(0) = Vq 



sinh" 2 (0/20) p < 0, 



(4.23) 



with Vq = 487r 3 a>/fi: 4 |/3|. 

Note that for 0/0 ^> 1, V(0) oc exp(— 0/0) for both choices of the sign of (3. 
Large values of cross the Hubble radius late in inflation and correspond to small A. 
Notice from Eq. (|4.19[) that Aq — > X v as A — >• 0. We have already reconstructed the 
potential that results from this A G {X) as V(0/0) oc exp(— 0/0), which agrees with 
the definition of d> given in Eq. ( gig ) when y<l. (The assumption that Aq As 
is equivalent to this condition). 
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We can also expand Eq. (|4.23|) for small 0: 



1 - (0/0) 2 /4 + • • • p>o 
V(<j>) — Vq { (4.24) 
4(0/0)- 2 + --- /3<0. 

The positive (3 case is also an approximation to a potential of the form V(0) oc 
1 + cos(0/0) as studied in a type of model called natural inflation . 



The purpose of the above reconstruction exercises is to demonstrate how the 
reconstruction process proceeds. We have reconstructed several popular inflationary 
potentials from knowledge of either the scalar or tensor perturbation spectrum. Before 
turning to the prospectus for actually determining As and A G from observational 
data, in the next section we discuss a "perturbative" approach in reconstruction of 
the potential. 



V. RECONSTRUCTING A PIECE OF THE POTENTIAL 



The reconstruction program described in the previous section is quite ambitious, 
as it depends upon knowledge of the functional forms of Aa(X) and/or As(A) over a 
range of A. In this section we will outline a less ambitious, but more realistic program. 
We will assume that we have information only about the scalar and tensor spectra 
(and their first and second derivatives) at a single scale Ao, and see what we can learn 
about the potential.0 This "perturbative" approach to reconstruction may be useful 
in the very near future [f2T|. 

If we know A G {Xo) and As(Xo) at some length scale Ao (which left the Hubble 
radius during inflation when the value of the scalar field was 0o), we can use Eqs. 
([H| and (|3TT3| ) to find that 



1 dX 

A d4> 

dA G (X) 



k A|(A ) - A 2 G {\ 



dX 



A=A 



A=A 



y/2 A S (X )A G (X ) 
Ao A|(A )-^(Ao)' 



(5.1) 



l/(0o) immediately follows from A G (Xo) and As(Xc 



V(<f> ) 



3A G (\o) - 



A%{X, 



o 



(5.2) 



10 In practice, observing the derivatives at a single point may be just as difficult as measuring the 
shape over a range of scales, though one might hope for adequate information to be obtained from 
a significantly smaller range of scales (and with more freedom to coarse-grain), perhaps even those 
accessible from a single experiment. 



21 



With the approximation that Ag(Xq) <C As(Ao), the expression simplifies to 



48tt 



^(Ao) 



l + O 



' AUK 
A|(A ) 



(5.3) 



Further reconstruction of the potential requires more than simply knowledge of 
the amplitudes of the scalar and tensor perturbations at Ao, we must know the first 
derivative, or the spectral index of the scalar spectrum at A = A : 



dAs(X) 



A S (X) dX 



1 — n 



A=A 



2A 



A=A 



1 - n 
2An 



(5.4) 



Using Eqs. ( |5.1| ) and ( |5.2| ) we can find the first derivative of the potential: 



VI 



dV{(f)) 



A=A 



167T 3 AUXq) 
V2k 3 A s (X ) 



7 - n - (5 - n 



A 2 G (X ) 
AUXo) 



(5.5) 



This expression also simplifies in the Ag(Xq) ^ ^s(Ao) limit: 



V'( 



16tt 3 A 3 G (X C 
V2k 3 A s (X 



-(7 -n ) 



l + O 



(5.6) 



Repeated differentiation of this expression with respect to <fi enables one to derive 
the higher derivatives. In principle, the potential can then be expanded as a Taylor 
series about the point 0o- The full expression for the second derivative is 



V"( 



Stt 3 AUX ) 
k 2 A 2 s (Xo) 
1 - n 



3A 2 G (X 



o 



7 



+L4|(A ) - AUXo)} 

a%(x ) 



7 - n - (5 - n ) 
-1 + 



n - (5 - n ) 



^(Ac 



A|(A ) 
A 2 S (X )J 



AUXo) 

[A 2 s(X ) - A G (Xo)] 



Xon'o 



-2(5 -n ) 



+(5-n )(l -n ) 



AUXo) AUXo) - AUXo) 
AUXo) 



AUXo) 



(5.7) 



where n' Q = dno/dXo- If one makes the approximation that Ag{Xq) <C As(Xq), it 
follows that this expression simplifies considerably: 



V' 



4vr 3 AUXo) 
AUXo) 

l + O 



4(n - 4) 2 ^(A ) - (1 - no) (7 - n )^ 5 (A ) 



x 



AUXo 
AUXo 



o 



drir, 



<im A 



(5.8) 



A=A / 
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Note that (1 — no) is in principle of the same order as A G j A\. 

It is hoped that observations will soon be of a sufficient standard to measure 
Ag(\q), As(Xo) and n at a particular point. One may then be able to establish 
whether the potential is convex or concave with the use of Eq. fl5.8|) . The potential is 
convex in any model where the scalar spectrum decreases with increasing wavelength. 
Note though that this is only sufficient, not necessary. Indeed, most popular models 
such as polynomial and exponential potentials are convex yet still feature a spectrum 
increasing with increasing wavelength. In models where the tensor contribution is 
negligible, the only important parameter is the sign of 1 — no, since no > 7 is already 
ruled out by observation 

Another quantity of interest that may soon be determined observationally is the 
dimensionful parameter V(<po) / \ V'((p ) \ . This is determined by the relative amplitudes 
of the scalar and tensor fluctuations at a given scale via Eq. ( |3.5|) . In this sense, such 
a quantity yields information regarding a mass scale at which these processes are 
occurring during inflation. In the case of polynomial potentials it uniquely determines 
(j)Q. For exponential and hyperbolic examples, however, it measures the steepness of 
the potentials as given by 0. 

Although the value of 0o is undetermined because of the inherent freedom to shift 
^ by a constant, some information of the range of covered by observations of the 
spectra on scales between Ao and Ai can be recovered. We can start with Eq. (|3.13| ), 
the reconstruction equation for 0(A) and find 



where 0i = 0(Ai). Then, we can then use a simple trapezoidal integration rule to 
find 0i — 0o in terms of the spectra at Aq and Ai. 



Conventionally, one chooses a particular theory, assesses the spectrum it predicts 
and attempts a comparison between its predictions and the observed Universe. For 
our purposes here, one must be more ambitious and execute this procedure in reverse 
order, proceeding from the observations to the primeval spectrum, and thence to 
the underlying inflationary theory. As well as covering the current observational 
position, we intend to survey the possibilities inherent in future experiments, both 
proposed and conjectural, in determining the primeval spectrum.0 In keeping with 
our inflationary motivation, we assume throughout that we have a universe of critical 
density. 

11 For an extensive review of large-scale structure studies, see the papers of Efstathiou j2j| and 
Liddle and Lyth §§. 




(5.9) 



VI. DETERMINING THE PRIMEVAL SPECTRUM 
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The range of scales of interest stretches from the present horizon scale, GOOO/i^ 1 
Mpc, down to about lh" 1 Mpc, the scale which contains roughly enough matter to 
form a typical galaxyQ On the microwave sky, an angle of 9 (for small enough 9) sam- 
ples linear scales of lOO/i -1 (#/l°)Mpc.0 For purposes of discussion, it is convenient 
to split this range into three separate regions. 

• Large scales: 6000/1" 1 Mpc — > ~ 200/i" 1 Mpc: 

These scales entered the horizon after the decoupling of the microwave back- 
ground. Except in models with peculiar matter contents, perturbations on these 
scales have not been affected by any physical processes, and the spectrum retains 
its original form. At present the perturbations are still very small, growing in 
the linear regime without mode coupling. Here, we are still seeing the primeval 
spectrum. 

• Intermediate scales: ~ 200/i -1 Mpc — ► 8/i 1 Mpc: 

These scales remain in the linear regime, and their gravitational growth is easily 
calculable. However, they have been seriously influenced by the matter content 
of the Universe, in a way normally specified by a transfer function, which mea- 
sures the decrease in the density contrast relative to the value it would have 
had if the primeval spectrum had been unaffected. Even in CDM models, where 
the only effect is the suppression of growth due to the Universe not being com- 
pletely matter dominated at the time of horizon entry, this effect is at the 25% 
level at 200/i _1 Mpc. To reconstruct the primeval spectrum on these scales, it 
is thus essential to have a strong understanding of the matter content of the 
Universe, including dark matter, and of its influence on the growth of density 
perturbations. 

• Small scales: 8/i 1 Mpc — ► Ihr 1 Mpc: 

On these scales the density contrast has reached the nonlinear regime, cou- 
pling together modes at different wavenumbers, and it is no longer easy to 
calculate the evolution of the density contrast. This can be attempted either 
by an approximation scheme such as the Zel'dovich approximation [E^IRO], or 



more practically via iV-body simulations |31|. Further, hydrodynamic effects 
associated with the nonlinear behavior can come into play, giving rise to an 
extremely complex problem with important non-gravitational effects. Again, 
the transfer function plays a crucial role on these scales. In hot dark matter 
models, perturbations on these scales are most likely almost completely erased 
by free streaming, and hence no information can be expected to be available 
(far less than the rather detailed information reconstruction would require). In 



12 The present density of the Universe is p = 3H$/8nG ~ 2.8k- 1 x 10 11 M Q (h~ 1 Mpc)~ 3 , where 
Mq is the solar mass. 

13 The surface of last scattering is located some 200/i -1 Mpc inside the horizon distance. 
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a CDM model, enough residual perturbations may remain on these scales for 
useful information to be obtained. 

Let us now consider each range of scales in turn, starting with the largest scales 
and working down to the smallest scales. 

A. Large scales (6000/r 1 Mpc — ► ~ 200/i" 1 Mpc): 

Without doubt the most important form of observation on large scales for the near 
future is large-angle microwave background anisotropics. Scales of a couple of degrees 
or more fall into our definition of large scales. Such measurements are of the purest 
form available — anisotropy experiments directly measure the gravitational potential 
at different parts of the sky, on scales where the spectrum retains its primeval form. 
Such measurements also are of interest in that the tensor modes may contribute. 
Tensor modes do not participate in structure formation and most measurements we 
shall discuss are oblivious to them. Further, tensor modes inside the horizon redshift 
away relative to matter, and so tensor modes also fail to participate in small-angle 
microwave background anisotropies. 

Nevertheless, these large-scale measurements still exhibit one crucial and ulti- 
mately uncircumventable problem. On the largest scales, the number of statistically 
independent sample measurements that can be made is small. Given that the un- 
derlying inflationary fluctuations are stochastic, one obtains only a limited set of 
realizations from the complete probability distribution function. Such a subset may 
insufficiently specify the underlying distribution, which is the quantity predicted by 
an inflationary model, for our purposes. This effect, which has come to be known 
as the cosmic variance, is an important matter of principle, being a source of un- 
certainty which remains even if perfectly accurate experiments could be carried out. 
At any stage in the history of the Universe, it is impossible to accurately specify 
the properties (most significantly the mean, which is what the spectrum specifies 
assuming gaussian statistics) of the probability distribution function pertaining to 
perturbations on scales close to that of the observable Universe. 

Observations other than microwave background anisotropies appear confined to 
the long term future. Even such an ambitious project as the Sloan Digital Sky Sur- 
vey (SDSS) |32| can only reach out to perhaps 500/i -1 Mpc, which can only touch 
the lower end of our specified large scales. However, in order to specify the fluctu- 
ations accurately, one needs many statistically independent regions (100 seems an 
optimistic lower estimate) which means that the SDSS may not specify the spectrum 
with sufficient accuracy above perhaps IOO/tT 1 Mpc. 

A much more crucial issue is that the SDSS will measure the galaxy distribution 
power spectrum, not the mass distribution power spectrum that is our inflationary 
prediction. In modern work it is taken almost completely for granted that these 
are not the same, and it seems likely too that a bias parameter (relating the two 
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by a multiplicative constant) which remains scale independent over a wide range of 
scales may be hopelessly unrealistic. Consequently, converting from the galaxy power 
spectrum back to that of the matter may require a detailed knowledge of the process 
of galaxy formation and the environmental factors around distant galaxies. Once 
one attempts to reach yet further galaxies with a long look-back time, one must also 
understand something about evolutionary effects on galaxies. As we shall discuss in 
the section on intermediate scales, it seems likely that peculiar velocity data may be 
rather more informative than the statistics of the galaxy distribution. 

A more useful tool for large scales is microwave background anisotropies on large 



angular scales. Our formalism closely follows that of Scaramella and Vittorio [33 |. 
On large angular scales, the most convenient tool for studying microwave background 
anisotropies is the expansion into spherical harmonics 

A rp OO / 

— (x, 9, 0) = £ £ a /m (x) Yi(6, 0), (6.1) 

1 1=2 m=-l 

where 9 and are the usual spherical polars and x is the observer position. With 



spherical harmonics defined as in [34], the reality condition is 



o,,_ m = (-ir^ |m . (6.2) 

In the expansion, the unobservable monopole term has been removed. The dipole 
term has also been completely subtracted; the intrinsic dipole on the sky cannot 
be separated from that induced by our peculiar velocity relative to the comoving 
frame, though it is easy to show that for adiabatic perturbations it will be negligible 
compared to it. 

With gaussian statistics for the density perturbations, the coefficients aj m (x) are 
gaussian distributed stochastic random variables of position, with zero mean and 
rotationally invariant variance depending only on I 

(a, m (x))=0 ; (h m (x)| 2 > = (6-3) 

It is crucial to note that a single observer such as ourselves sees a single realization 
from the probability distribution for the ai m . The observed multipoles as measured 
from a single point are defined as 

Q 2 i=j- E M 2 > ( 6 - 4 ) 

m=—l 

and indeed the temperature autocorrelation function can be written in terms of these 

/AT AT \ 00 

C(a) = (— {9 1 ^ l )—{9 2 , fa)) a = E^(cosa), (6.5) 
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where the average is over all directions on a single observer sky separated by an angle 
a, and Pj(cosa) is a Legendre polynomial. The expectation for the Q 2 , averaged over 
all observer positions, is just 

(Q 2 ) = i-(2Z + 1)E 2 . (6.6) 

A given model predicts values for the averaged quantities (Q 2 )- On large angular 
scales, corresponding to the lowest harmonics, only the Sachs- Wolfe effect operates. 
One has two terms corresponding to the scalar and tensor modes — we denote these 
contributions by square brackets. The scalar term is given by the integral 

E?[5] = %T f jf (2k/aE) - 2 Amk) } (6.7) 
mr Jo k m z 

where ji is a spherical Bessel function and the transfer function T(k) is normalized to 
one on large scales. As an example, a power-law spectrum kA 2 s oc k n on scales where 
the transfer function is sufficiently close to unity gives the oft-quoted 

y2[91 _ y2[91 r(i + (n-l)/2) r((9-n)/2) 

^ [bl ~ ^ r(Z + (5-n)/2) r((3 + n)/2)' [ } 

which for a flat n — 1 spectrum gives the simple E 2 [S] = 6E 2 .[S']/7(Z + 1). However, 
true reconstruction requires the integral expression. 

The equivalent expression for the tensor modes is a rather complicated multi- 
ple integral which usually must be calculated numerically [fT4] , |15| ,pO|l . Under many 
circumstances (Lucchin, Matarrese and Mollerach |20| suggest 0.5 < n < 1 for power- 
law inflation) there is a helpful approximation which is that the ratio [5]/E^[T] is 
independent of I and given by 

m = ±^k (69) 

For many purposes this is a perfectly adequate expression, but for true reconstruction 
of the inflaton potential, one must of course use the exact integral expression. 

On the sky, one does not observe each contribution to the multipoles separately. 
As uncorrelated stochastic variables, the expectations add in quadrature to give 

S 2 = S 2 [5] + S 2 [T]. (6.10) 

For reconstruction purposes, there are two obstructions of principle. These are 

• Even if one could measure the E 2 exactly, the last scattering surface being closed 
means one obtains only a discrete set of information — a finite number of the 
Hi covering some effective range of scalesQ There will thus be an uncountably 
infinite set of possible spectra which predict exactly the same set of Ej. 



14 The l-th multipole is often taken as corresponding roughly to a scale ki ~ IHq /2Mpc 1 = 
Z/i/SOOOMpc- 1 . 
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• One cannot measure the Hf exactly. What one can measure is a single realiza- 
tion, the Qf. As a sum of 21 + 1 gaussian random variables, Qf has a probability 
distribution which is a x 2 distribution with 21 + 1 degrees of freedom, X21+V The 
variance of this distribution is given by 

Var[Q?] = i^(Q 2 }\ (6.H) 

though one should remember that the distribution is not symmetric. Each Qf 
is a single realization from that distribution, when we really want to know the 
mean. From a single observer point, there is no way of obtaining that mean, and 
one can only draw statistical conclusions based on what can be measured. Thus, 
a larger set of spectra which give different sets of T,f can still give statistically 
indistinguishable sets of Qf. The variance falls with increasing I but is significant 
right across the range of large scales. This is illustrated in Figure 2. 

Finally, it should be mentioned that measurements of the polarization of the 
CMBR on large scales may allow a separate determination of the gravitational wave 



spectrum to be made Such an effect is potentially detectable if gravitational 

waves dominate the COBE result and the polarization is of the order of 10%, say. 
If the waves only contribute 10% of the COBE signal, for example, then only 10% 
of 10% is polarized, which significantly reduces the overall effect. Unfortunately, 
reconstruction of the potential must await a positive detection of such an effect, so 
we will not discuss it further. 



B. Intermediate scales (~ 200/i 1 Mpc — > 8/T 1 Mpc): 



It is on intermediate scales that determination of the primeval spectrum is most 
promising, though sadly these scales only encompass about 3 e-foldings. Here a range 
of promising observations are available, particularly towards the small end of the range 
of scales. In terms of technical difficulties in interpreting measurements, a trade-off 
has been made compared to large scales; on the plus side, the cosmic variance is a 
much less important player as far more independent samples are available, while on 
the minus side the spectrum has been severely affected by physical processes and thus 
has moved a step away from its primeval form.0 



1. Intermediate-scale microwave background anisotropies 

15 In the distant future, when the horizon size is vastly greater than at the present, there will be 
a range of scales above 200ft. -1 Mpc where the cosmic variance remains small and the spectrum 
retains it primeval form. Such a region would be an ideal place to carry out reconstruction, but 
unfortunately does not exist at the present epoch. 
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In the absence of reionization, the relevant angular scales are from about 2° down 
to about 5 arcminutes. (Should reionization occur, a lot of the information on these 
scales could be erased or amended in difficult to calculate ways.) Several experiments 
are active in this range, including the South Pole and MAX experiments, but as none 
have yet published a positive detection they are not of direct interest to us here at 
present. Indeed, even with a detection many of these ground based experiments are 
unable to give results with the statistical quality we would require due to the small 
sky coverage which is typically involved. 

Unlike the large-scale anisotropy, one cannot write down a simple expression for 
the intermediate-scale anisotropics, even if it is assumed that one has already incorpo- 
rated the effect of dark matter on the growth of perturbations via a transfer function. 
The reason is due to the complexity of physical processes operating. A case in point 
is the expected anisotropy (specified by the £ 2 , but now for larger I) in the CDM 



model (n = 1), as calculated in detail by Bond and Efstathiou ^6 . 

On large scales, Z 2 E 2 is approximately independent of I as we have seen. Once 
we get into the intermediate regime, / 2 £ 2 exhibits a much more complicated form, 
which is dominated by a strong peak at around / = 200. This is induced by Thomson 
scattering from moving electrons at the time of recombination. Bond and Efstathiou's 
calculation gives a peak height around 6 times as high as the extrapolated Sachs- Wolfe 
effect. Beyond the first peak is a smaller subsidiary peak at / ~ 800. 

In their calculation, Bond and Efstathiou assumed both the primeval spectrum 
and the form of the dark matter. For reconstruction purposes, it seems that a good 
knowledge of the form of dark matter is a pre-requisite, in order that these processes 
can be calculated at all. Of course, given the number of active and proposed dark 
matter search experiments, one should be optimistic that this information will be 
obtained in the not too distant future. However, even with this information, the 
complexity of the calculation makes it hard to conceive of a way of inverting it, 
should a good experimental knowledge of the S 2 (I G [30,750]) be obtained. Once 
again, it's much easier to compare a given theory with observation than to extract a 
theory from observation. 

One of the interesting applications of these results might be in combination with 
the large-scale measurements. The peak on intermediate scales is due only to pro- 
cesses affecting the scalar modes, whereas we have pointed out that the large-scale 
Sachs-Wolfe effect is a combination of scalar and tensor modes. On large scales, 
one cannot immediately discover the relative normalizations of the two contributions. 
However, if the dark matter is sufficiently well understood, the height of the peak 
in the intermediate regime gives this information. Should it prove that the tensors 
do play a significant role, then this would be a very interesting result as it immedi- 
ately excludes slow-roll potentials for the regime corresponding to the largest scales. 
Should the tensors prove negligible, then although the conclusions are less dramatic 
one has an easier inversion problem on large scales as one can concentrate solely on 
scalar modes. 
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2. Galaxy clustering in the optical and infrared 



A. Redshift surveys in the optical. 

Over the last decade, enormous leaps have been made in our understanding of the 
distribution of galaxies in the Universe from various redshift surveys. Most prominent 
is doubtless the ongoing Center for Astrophysics (CfA) survey |37[| , which aims to form 
a complete catalogue of galaxy redshifts out to around 100/i -1 Mpc. Other surveys 
of optical galaxies, often trading incompleteness for greater survey depth, are also in 
progress. On the horizon is the Sloan Digital Sky Survey |32] which aims to find the 



redshifts of one million galaxies, occupying one quarter of the sky, with an overall 
depth of 500/i _1 Mpc and completeness out to 100/i -1 Mpc. 

The redshifts of galaxies are relatively easy (though time consuming) to measure 
and interpret, and so provide one of the more observationally simple means of de- 
termining the distribution of matter in the Universe. The main technical problem 
is to correct the distribution for redshift distortions (which gives rise to the famous 
"fingers of God" effect). However, the distribution of galaxies, specified by the galaxy 
power spectrum (or correlation function) is two steps away from telling us about the 
primeval mass spectrum. 

• We have already discussed that the primeval spectrum on intermediate scales 
has been distorted by a combination of matter dynamics and amendments to 
the perturbation growth rate when the Universe is not completely matter dom- 
inated. If we know what the dark matter is, then this need not be a serious 
problem. 

• Galaxies need not trace mass, and in modern cosmology it is almost always 
assumed they do not. This makes the process of getting from the galaxy power 
spectrum to the mass power spectrum extremely non-trivial. Models such as 
biased CDM rely on the notion of a scale-independent ratio between the two, but 
this too can only be an approximation to reality. In recent work, authors have 
emphasised the possible influence of environmental effects on galaxy formation 
(for instance, a nearby quasar might inhibit galaxy formation H38|1 ), and indeed 
it has been demonstrated that only very modest effects are required in order to 
profoundly affect the shapes of measured quantities such as the galaxy angular 
correlation function [091. 



Despite this, attempts have been made to reconstruct the power spectrum from 
various surveys. In particular, this has been done for the CfA survey |4(|, and for 



the Southern Sky Redshift Survey These reconstructions remain very noisy, 

especially at both large scales (poor sampling) and small scales (shot noise and redshift 
distortions), and at present the best one could do would be to try and fit simple 
functional forms such as power-laws or parametrized power spectra to them. Even 
then, the constraints one would get on the slope of say a tilted CDM spectrum are 
very weak indeed. However, these reconstructions go along with the usual claim that 
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standard CDM is excluded at high confidence due to inadequate large-scale clustering, 
without providing any particular constraints on the choice of methods of resolving 
this conflict. 

Nevertheless, with larger sampling volumes such as those which the SDSS will 
possess, one should be able to get a good determination of the galaxy power spectrum 
across a reasonable range of scales, perhaps 10/i -1 to 100/i -1 Mpc. 

B. Redshift surveys in the infrared. 

A rival to redshifts of optical galaxies is those of infra-red galaxies, based on 
galaxy positions catalogued by the Infra-Red Astronomical Satellite (IRAS) project 
in the mid-eighties. The aim here is to sparse-sample these galaxies and redshift the 
subset. This is being done by two groups, giving rise to the QDOT survey f|2]| and 



the 1.2 Jansky survey ||43| . Taking advantage of the pre-existing data-base of galaxy 
positions has allowed these surveys to achieve great depth with even sampling and 
reach some interesting conclusions. 

The main obstacle to comparison with optical surveys is due to the selection 
method. Infra-red galaxies are generally young, and appear to possess a distribu- 
tion notably less clustered than their optically selected counterparts. They are thus 
usually attributed their own bias parameter which differs from the optical bias. The 
mechanics of proceeding to the power spectrum are basically the same as for optical 
galaxies. 

The most interesting and relevant results here are obtained in combination with 
peculiar velocity information, as discussed below. 

C. Projected catalogues. 

As well as redshift surveys, one also has surveys which plot the positions of galax- 
ies on the celestial sphere. At present the most dramatic is the APM survey [03 



encompassing several million galaxies. The measured quantity is the projected coun- 
terpart of the correlation function, the angular correlation function usually denoted 
w(8) where 9 is the angular separation. Though arguments remain as to the presence 
of systematics, one in principle has accurate determinations of the galaxy angular cor- 
relation function. The first aim is to reconstruct the full three dimensional correlation 
function from this (proceeding thence to the galaxy power spectrum). Unfortunately, 
present methods of carrying out this inversion (based on inverting Limber's equa- 
tion which gives w{9) from £(r) ) have proven to be very unstable, and a satisfactory 
recovery of the full correlation function has not been achieved. 

In its preliminary galaxy identification stage, the SDSS will provide a huge pro- 
jected catalogue on which further work can be carried out. 

3. Peculiar velocity flows 
Potentially the most important measurements in large-scale structure are those 
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of the peculiar velocity field. Because all matter participates gravitationally, peculiar 
velocities directly sample the mass spectrum, not the galaxy spectrum. Were one to 
know the peculiar velocity field, this information is therefore as close to the primeval 
spectrum as is microwave background information. Indeed, in the linear regime the 
spectrum of the modulus v of the velocity^ is just given by 

Perhaps the most exciting recent development in peculiar velocity observations is 
the development of the POTENT method by Bertschinger, Dekel and collaborators 
45fl . Using only the assumption that the velocity can be written as the divergence 



of a scalar (in gravitational instability theories in the linear regime this is naturally 
associated with the peculiar gravitational potential), they demonstrate that the radial 
velocity towards/away from our galaxy (which is all that can be measured by the 
methods available) can be used to reconstruct the scalar, which can then be used to 
obtain the full three dimensional velocity field. This has been shown to work very 
well in simulated data sets, where one mimics observations and then can compare 
the reconstruction from those measurements with the original data set. So far, the 
noisiness and sparseness of available real radial velocity data has meant that attempts 
to reconstruct the fields in the neighborhood of our galaxy have not yet met with great 
success; however, once better and more extensive observational data are obtained one 
can expect this method to yield excellent results. 

At present, POTENT appears at its most powerful in combination with a sub- 
stantial redshift survey such as the IRAS/QDOT survey. As POTENT supplies 
information as to the density field and the redshift survey to the galaxy distribution, 
the two in combination can be used in an attempt to measure quantities such as the 
bias parameter and the density parameter Q of the Universe. Reconstructions of 
the power spectrum have also been attempted |[46|| . At present, the error bars (due 
to cosmic variance because of small sampling volume, due to the sparseness of the 
data in some regions of the sky and due to iterative instabilities) are large enough 
that a broad range of spectra (including standard CDM) are compatible with the 
reconstructed present-day spectrum. 

With larger data sets and technical developments in the theoretical analysis tools, 
POTENT (and indeed velocity data in general) appears to be a very powerful means 
of investigating the present-day power spectrum. To that, one need only add a knowl- 
edge of the dark matter to obtain the primeval spectrum and thence to the inflaton 
potential. Although likely to be limited to the range of scales specified at the lower 
end by the onset of the nonlinear regime and at the upper end by the range of feasi- 
ble experimental measurements of the radial peculiar velocity, it seems that velocity 



16 The spectrum is defined as V v — V(k 3 /2ir 2 )(\5 v \ 2 } , with V being the volume over which the 
Fourier components 5 v (k) are defined. 
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data provide the most promising means of reconstructing a segment of the inflaton 
potential. 



C. Small scales (8/i 1 Mpc — ► Ik' 1 Mpc): 

It is worth saying immediately that this promises to be the least useful range of 
scales. For many choices of dark matter, including the standard hot dark matter 
scenario, perturbations on these scales are almost completely erased by dark matter 
free-streaming to leave no information as to the primeval spectrum. Only if the dark 
matter is cold does it seem likely that any useful information can be obtained. 

There are several types of measurement which can be made. Quite a bit is known 
about galaxy clustering on small scales, such as the two-point galaxy correlation 
function. However, the strong nonlinearity of the density distribution on these scales 
erases information about the original linear-regime structure, and the requirement 
of TV-body simulations to make theoretical predictions makes this an unpromising 
avenue for reconstruction even should nature have chosen to leave significant spectral 
power on these scales. There exist very small-scale (arcsecond-arcminute) microwave 



background anisotropy measurements | 47| , though these are susceptible to a number 



of line of sight effects, and further the anisotropics are suppressed (exponentially) 
on short scales because the finite thickness (about 7h~ l Mpc) of the last scattering 
surface comes into play. 

Up to now, the most useful constraints on small scales have come from the pair- 
wise velocity dispersion |IS| (the dispersion of line-of-sight velocities between galax- 
ies). These are sensitive to the normalization of the spectrum at small scales, though 
unfortunately susceptible to power feeding down from higher scales as well. There are 
certainly noteworthy constraints — for instance it is generally accepted that unbiassed 
standard CDM generates excessively large dispersions. However, the calculations 
required involve iV-body simulations and because a wide range of wavelengths con- 
tribute, obtaining knowledge of any structure in the power spectrum on these scales 
is likely to prove impossible, even if the amplitude can be determined to reasonable 
accuracy. 



VII. DISCUSSION AND CONCLUSIONS 

To date, the traditional approach in cosmology has been to take a set of theoretical 
predictions for the structure of our universe and compare them directly with what is 
observed. The aim is to reduce to a minimum the space of possible theories consistent 
with observations. Unfortunately such an analysis can only deduce which theories are 
unsuitable and is unlikely to select uniquely the correct one. An alternative and more 
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ambitious program is to use the observations to reconstruct the theory. Within the 
context of the inflationary universe, for example, such an approach is justified when 
one considers the prize on offer — the form of the inflaton potential. The purpose of 
the present work has been to illustrate how such a reconstruction of the potential is 
possible in principle. 

There are two steps to any reconstruction procedure. In practice the observational 
information may not be in a form which allows a direct comparison with the theoret- 
ical predictions to be made. It is therefore necessary to first convert the data into the 
quantity predicted and only then can the second step of reconstructing the potential 
be completed. As was shown in Section VI, this is especially true in the inflationary 
universe and presents a number of fundamental difficulties with the procedure. 

In Section III, however, we successfully completed the second step of the process 
by deriving the correspondences between the tensor and scalar fluctuation spectra and 
the potential. This extended the analysis of Ref. [0] and a number of examples were 
presented in Section IV. In a true reconstruction one should make no assumptions 
concerning the form of the potential. In particular, the assumptions of slow-roll, which 
are essentially conditions on the flatness of V(<p), should be avoided. The formalism 
used places no restrictions on the inflaton field dynamics, but does assume the slow- 
roll expressions for the perturbation spectra still apply. From a computational point 
of view, it follows that reconstruction is unambiguous once the tensor spectrum is 
known. Unfortunately, however, it is this quantity which is the most difficult to 
determine observationally. The only observational effect of primordial gravitational 
waves appears to be their influence on large-scale CMBR anisotropics. We conclude 
that the most promising method of determining the tensor spectrum is to combine the 
large-scale CMBR results with intermediate scale data from peculiar velocities and 
CMBR anisotropies. The latter require a knowledge of the dark matter component, 
but are independent of any bias in the galaxy distribution. They determine the scalar 
spectrum, whereas the former depends on both the scalar and tensor modes. A simple 
subtraction therefore yields the tensor spectrum. 

Eq. ( p. 11 ) will allow a test of the inflationary paradigm to be made if a separate 
determination of the tensor spectrum on large scales can be made. A separate deter- 
mination of Aq on large angular scales coupled with COBE [|I|], Tenerife [^] and the 
Princeton-MIT balloon |5(J would lead to A 2 S . This could then be compared with the 
theoretical prediction derived from Eq. ( |3.11| ). If a discrepancy was found, it would 
suggest that one or more of the initial assumptions — such as the background space- 
time being flat; using a single, minimally coupled scalar field or Einstein gravity — were 
incorrect. On the other hand, in the absence of any discrepancy, this result could be 
used with a combination of CMBR measurements around 2°, velocity and galaxy clus- 
tering data, and compared with the theoretical predictions for different dark matter 
models. This would lead to limits on the form of dark matter present in the Universe. 

We note that reconstruction is still possible if the gravitational waves are not 
significant, although one must then deal with the integration constant which arises 
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in the solution of Eq. ( |3.11|) and can affect the functional form of the potential. 

Although we have been somewhat pessimistic about the near-term prospects for 
reconstructing the functional form of the potential, we are optimistic regarding the 
near-term possibility of obtaining some knowledge about the potential. To illustrate 
the promise of our method, let us assume that within a few years that a combination 
of CMBR measurements give us some information about the scalar and tensor am- 
plitudes at a particular length scale A (corresponding to an angular scale 9 ). An 
example is that we might in the near future have in hand the following: 

A S (X ) = 1 x 10- 5 ; A G (X ) = 2 x 1(T 6 ; n = 0.9; n' = 0. (7.1) 

If we would have this information, we can follow the perturbative procedure out- 
lined in Section V and reconstruct information about the potential in the vicinity of 
some point 4>o'- 

V{<f) ) = (2 x 10 16 GeV) 4 
±V'(4> Q ) = (3 x 10 15 GeV) 3 
V"((f) ) = (5 x 10 13 GeV) 2 . (7.2) 

By taking some appropriate ratios one may find mass scales for the potential. In this 
way cosmology might be first to get a "piece of the action" of GUT-scale physics. 
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Figure Captions 

Figure 1: 

A schematic figure illustrating the main concepts behind reconstruction. For infla- 
tion the two main steps involve converting the observations (lower half of figure) into 
the primordial scalar (As) and tensor (Aq) fluctuation spectra and then working in 
reverse to reconstruct the potential V ((/)). The main observational information from 
the cosmic microwave background arises through the Cosmic Background Explorer 
(COBE) satellite 0, and the Tenerife (TEN) g9[ and South Pole (SP) f5TJ col- 
laborations. Galaxy surveys (APM g|, CfA |37J, IRAS may provide useful 
information up to lOO/iT 1 Mpc, while the Sloan Digital Sky Survey (SDSS) |52) should 
extend to the lowest scales measured by COBE. Peculiar velocity measurements using 
the POTENT (P) |^5[ methods are important on intermediate scales. The angle 9 
measures angular scales on the CMBR in degrees, and length scales A are in units of 
h^ 1 Mpc. dn refers to the horizon size today and at recombination and g?nl ~ 8/i _1 
Mpc is the scale of non-linearity. (See the text for details). Perfect observations will 
only reconstruct a small portion of the inflaton potential corresponding to between 
53 < AN < 60 e-foldings before the end of inflation. 



Figure 2: 

Multipoles up to Z = 30, roughly corresponding to the complete range of large scales. 
The solid line represents the ensemble averaged (Qf) (multiplied by I) for a flat 
(n = 1) spectrum of scalar density perturbations with A G (X) <C A S (X) , normalized 
to Eg = 1. The three dashed lines represent different randomly chosen realizations 
of this distribution. Observations can only supply one such line, giving little clue 
to the ensemble average quantity that inflation supplies the form of. For compari- 
son, the dash-dotted line shows the result of a scalar spectrum with n = 0.8, again 
with Ag{\) <C As (A) (such a combination would arise from an appropriate inverted 
harmonic oscillator potential). Note that the normalization of this line is arbitrary 
(shown here with E| = 1), and were it moved up it could match an observed distribu- 
tion across much of the range. More significantly, it is easy to note that any detailed 
information in the spectrum such as peaks or troughs can be swamped completely by 
the cosmic variance. 
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